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Asteroid Threat Assessment Project 


Planetary Defense Coordination Office 
Mission Statement: 


This new office was recently established at NASA HQ to coordinate 
planetary defense related activities across NASA. 


Lead national and international efforts to: 


¢ Detect any potential for significant impact of planet Earth by 
natural objects 


¢ Appraise the range of potential effects by any possible impact 
¢ Develop strategies to mitigate impact effects on human welfare 
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ATAP’s Overarching Assignment 


¢ Provide information needed to determine the minimum size asteroid 
for which in space mitigation must be undertaken, and therefore must 
be detected in time to take action. 


¢ Provide information needed to determine the maximum size asteroid 
for which civil defense measures are sufficient for mitigation actions. 
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2013 Chelyabinsk Event 


In February 2013, the blast over Chelyabinsk in Russia released the equivalent of 
500,000 tonnes of TNT. Around 1,000 people were injured in that explosion - 
mostly as result of flying glass from smashed windows. The meteoroid, estimated 
to be about 10 tons, entered the atmosphere at a speed of 19.15 km/s. 
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Introduction 


¢ Comet - Cosmic snowballs of frozen gases, rock 
and dust roughly the size of a small town. 


¢ Asteroid - A small rocky body orbiting the sun. 
Large numbers of these, ranging in size from 
nearly 1,000 km to dust particles, are found as the 
asteroid belt. 


¢ Meteoroid - a small rocky or metallic body in 
outer space that range in size from small grains to 
1 meter-wide objects. 


¢ Bolides - an extremely bright asteroid, that 
explodes in the atmosphere. 


¢ Meteorites - A meteorite is a solid piece of debris 
from an object, such as a comet, asteroid, or 
meteoroid, that originates in outer space and 
Survives its passage through the Earth's 
atmosphere and impact with the Earth's surface. 


Asteroid 


Bolide 
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Background 


¢ The asteroids can be classified into following broad categories 
— C-type (Carbonaceous) 
— $-type (Stony) 
— X-type (Mostly Metallic) 
¢ The size of an entering asteroid could range from a few meters to 
several kilometers 
¢ Very high entry speed leading to high stagnation pressures, presence 
of pre-existing cracks cause some of the asteroids (esp. in the S-type 
category) to break during entry. 
¢ The fragmentation could be a chain event and could occur at multiple 
altitude. 
¢ For entry break-up and modeling activities we first focused on stony 
asteroids in less than 100m range. 
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Entry Capsules vs. Asteroid Entry 


Apollo Entry Chelyabinsk Entry 


¢ Can some of the modern computational analysis tools used in 
design of Entry capsules be used for simulation of asteroid 
entries? 
¢ Can we develop asteroid models for : 
— Material thermal response 
— Material structural response, including fragmentation 


- Energy deposition along asteroid trajectory in the atmosphere, i.e., 
light curves 
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Entry Systems Tool >Asteroid Entry 


Shape Regular and smooth geometry Irregular and rough geometry 
Highly complex mineralogy - 
Manufactured ablative material Depends on asteroid class (S, M, 


Material 
X) 
~ . Could have voids & could be 
Minimal voids and cracks; ; 
Structure fractured; 


SCTE SMe 2 Welt) ee nteks Structural properties of meteorites 


Shape Change | Negligible recession/mass loss Recession dominated 
Fragmentation — | Not an option! S-class most likely to fragment 


Capsule entry physics has some things in common with meteor 
physics, but approaches to the problem are different — prediction vs 


reconstruction 
Bringing reliable predictive capabilities to bear on asteroid entries is the 


focus of the effort. 
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Fragmentation Modeling Objectives 


To understand fragmentation and fracture of a given asteroid and 
mechanisms of break-up. 
¢ Determine mechanical properties of asteroids 
¢ Develop finite element modeling techniques for stony asteroids 
¢ Investigate changes in the stress distribution in the presence of pre- 
existing defects 
¢ Predict fragmentation and release of energy as an asteroid is entering 
into Earth’s atmosphere for various sizes 


— 
PS. : OCTOBER 23 — 27, 2016 
ase eer SALT PALACE CONVENTION CENTER | SALT LAKE CITY, UTAH USA 


te - | 


Lay 


Methods of Validation 


¢ Ground testing of meteorites 


— These represent the strongest part of an incoming asteroid and pose 
challenge of scaling to asteroids 


¢ Ground observations 
— Light curves 
— Infrasound 
— Observed falls 


Challenge : Establish a link between meteorites and observations of 
asteroid entry 


Could we learn from terrestrial rocks and geophysics discipline? 
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Common form of large 
scale rock structures (rock 
masses) behave differently 
from small scale rock 
samples. 


There are features (cracks, 
flaws) that create some 
“average, macroscopic 
strength” 


The strength is size 
dependent 


Account for size” by 
natetettadlaleMdiccvalenia 
(Credit: K Holsapple) 
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Meteorites to Asteroids 


We plan to establish a two prong approach to establish the link 
between meteorite properties and full scale modeling of asteroids 


1. Develop meteorite material models (mechanical properties, constitutive 
relations and strength models ). Feed these models to full scale 
simulations of Asteroids 


2. Develop Weibull parameters based on ground tests, CT scans and FE 
simulations of meteorite models. Predict strength/corresponding 
fragmentation altitude of incoming Asteroid based on Weibull 
parameter. 
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Asteroid Entry Simulations — ALE3D 


Time: 3.842 s Altitude: 36.98 km 
Temperature (K) 
— 20000. 
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Material models — modulus, density, strength etc. are required for a 
high fidelity simulation of asteroid entry 


D. Robertson 


—- wa 
S OCTOBER 23 — 27, 2016 
SSS SALT PALACE CONVENTION CENTER | SALT LAKE CITY, UTAH USA 


> 7 


Scaling: Meteorites to Asteroids 


Bonding strength: 
Ground Truth ~10 cm Cracks, Veins 
Meteorites 


Terrestrial Analogue 


Account for size, bys 
modifying strength? 
<(Credit:.K Holsapple) 
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voids, grain size) Ex angle and breakup i | one might select off the ground for testing. 

3. Fusion crust (ablation) § | altitudes known. Ties | 

4. Mechanical properties = (to apparent 
(density, modulus, = “strength” 
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Development of Weibull Parameters 


; “modeling of full 
scale asteroid 


1. 10+ #100 ©1000 ~=10000 
Di(cm) 


Background figure taken from : Scale-dependent measurement of meteorite strengths (Desiree et. al.) 
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Properties of Meteorites 


¢ Most meteorites are very complex composites of following constituents 


— Large grains of several metallic and non-metallic minerals in form of 
chondrules 


— Porosity 
— Loose fine grains — matrix 
— Voids, cracks, metallic phase filled veins 
¢ Extensive literature survey was conducted to investigate material 
properties (Young’s modulus, Poison's ratio etc., fracture models) for 
stony meteorites. 
¢ However, we didn’t find any significant data on meteoritic properties. 
¢ [It became clear to us that in order to account for complex mineralogy 
and morphology of meteorites we will have to develop our own material 
models. 
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Meteorites vs. Terrestrial Analogs 


¢ The next step was to investigate terrestrial rock analogs like 
Basalt, Quartz and Gabbro. 
¢ Further investigations with characterization team revealed 


— The meteorites have several more minerals (metallic ones) that are not 
present in terrestrial rocks. 

— Presence of the metallic minerals significantly influences the 
properties. 


— In general the petrology (mineralogy as well as morphology) of stony 
meteorites are very different from terrestrial analogs 


© geology.com 


Meteorite 


Basalt 
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Meteorite Unit Model 


¢ Meteorite units are developed by randomly distributing the various 
minerals as different constituents in a 1 cm cube. 

¢ The Meteorite unit is divided into several small cubes representing the 
different constituents (smaller than chondrules) 

¢ Monte Carlo simulations are performed on several 100s of these units 
to distribute the constituents into many random orientation 

¢ These meteorite units represent the effective material properties taking 
into account the different compositions as well as petrology of the 
ordinary chondrites. 

¢ These effective properties can be used for both FE simulations as well 
as particle based models such as ALE-3D 
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Meteorite Unit Model - Approach 


¢ Three different meteorite categories are identified for meteorite unit 
development: 
— L (Low Iron — Ordinary Chondrite) 
— LL (Low Iron, Low metal Ordinary Chondrite) 
— H (high iron Ordinary Chondrite) 
¢ The unit models are also developed for rock analogs like basalt that 
can serve as a means to validate these models. 
¢ CT scans will be performed on ground based meteorites to 
investigate occurrence of micro and macro-cracks and other 
features like voids. 
— The results from these scan will determine whether and how to take these 
features into account when developing the constitutive relations. 
¢ Uniaxial compression tests on stony meteorites and rock analogs 
like basalt are performed to validate the meteorite unit models. 
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Mineral Composition — Stony Meteorites 


Principal minerals in 
chondritic meteorites (wt%) 


H oe a L (Ordinary 


Olivine Se ee ee 
Forsterite Mg2sioa =| 3.275 || | 8.245 | 85.156 | 
Fayalite fFe2sioa = || S439 | S| S55 | t.844 

Pyroxene (Mayreysiog _{___|_so0_|__a3-27 _|_24-25-_ 

Enstatite Mgsios—— (<s|| | —s—idL:CsitSs|sst7.779 CY 

Ferrosilite |FeSio3_ | = 3.95 | S| | 807 
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Different Phases — Stony Meteorites 


Table 1 Abundances of refractory inclusions, chondrules, metallic Fe,Ni, and matrix and other key properties of the chondrite groups. 


Group Refract. lith./Mg CAI and AOA Chondrule Chondrules Metal Matrix Fall frequency“ Examples 
rel. CI“ (vol. %) average diameter (vol.%)? (vol.%) (vol. %)° (%) 
(mm) 

Carbonaceous 

cI 1.00 <0.01 none <5 <0.01 95 0.5 Ivuna, Orgueil 

CM 1.15 5 0.3 20 0.1 70 1.6 Murchison 

co 1.13 13 0.15 40 1-5 30 0.5 Ornans 

CV 1.35 10 1.0 45 0-5 40 0.6 Vigarano, 
Allende 

CR 1.03 0.5 0.7 50-60 5-8 30-50 0.3 Renazzo 

CH 1.00 0.1 0.02-0.09 ~70 20 5 0 ALH 85085 

CB, 1.0 <0.1 ~5 40 60 <5 0 Bencubbin 

CB, 14 <0.1 ~0.5 30 70 <5 0 QUE 94411 

CK 1.21 4 0.8 15 <0.01 tie 0.2 Karoonda 

Ordinary 

H 0.93 0.01-0.2 0.3 60-80 8 10-15 34.4 Dhajala 

i 0.94 <0.1 0.5 60-80 3 10-15 38.1 Khohar 

LL 0.90 <0.1 0.6 60-80 1.5 10-15 7.8 Semarkona 

Enstatite 

EH 0.87 <0.1 0.2 60-80 8 <0.1-10 0.9 Qingzhen, Abee 

EL 0.83 <0.1 0.6 60-80 15 <0.1-10 0.8 Hvittis 

Other 

K 0.9 <0.1 0.6 20-30 6-9 70 0.1 Kakangari 

R 0.95 <0.1 04 >40 <0.1 35 0.1 Rumuruti 


Sources: Scott et al. (1996); other data from Weisberg et al. (1996, 2001), Rubin (2000), Krot et al (2002a), Kimura er al. (2002), and Bischoff et al (1993). 
* Mean ratio of refractory lithophiles relative to magnesium, normalized to Cl chondrites. * Includes chondrule fragments and silicates inferred to be fragments of chondrules. “Includes matrix-rich clasts, which account 
for all matrix in CH and CB, chondrites (Greshake et al., 2002). 4 Fall frequencies based on 918 falls of differentiated meteorites and classified chondrites (Grady, 2000). 
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H-type : Chemical constituents by volume 


The wt% were converted to volume % 


Mineral Model : 15% Matrix Volume Mineral Model: 10% Matrix Volume 


= Olivine 


m= Pyroxene 


NS = Diopside 
in = Feldspar 


= Troilite 
= Kamacite 
m@ Taenite 
m= Matrix 
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= Olivine 
m Pyroxene 
= Diopside 
= Feldspar 
= Troilite 
= Kamacite 
m Taenite 
m Matrix 


Meteorite Model — Material Map 
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Mechanical Properties of Meteorite Units 


¢ A 3.0 Mpa (30 atm) pressure 
applied based on trajectory 
analysis 

¢ The average strain tensor, e, 
provides the Young's modulus 
and Poisson's ratios of the 


cube. 
E, = O,/E,, 
Voy = Ey) Ez, 
y= E/E O,, Distribution in a cell 
ZX X. ZZ 
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Parameters Affecting the Unit Cell Predictions 


¢ Modulus of Olivine and Pyroxene 
¢ Porosity 
¢ % volume of Matrix material 
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Mineral Mechanical Properties 


Old (Higher) value of Young’s Modulus?* New (lower) values of Young’s Modulus? 


400 
350 Reference #1 ey Reference #2 
5 Olivine os 
© oe e 300 = Olivine 
& meee 9 m™ Pyroxene 
3 oa = Diopside 4 250 gos 
2 2 = Diopside 
3 200 m Feldspar 3 20 si 
eldspar 
= 150 = Triolite = i. Sp 
eo wi ; 
5 m Kamacite bp = Triolite 
(eo) Ss A 
> 100 m Taenite 2 100 = Kamacite 
50 a mw MATRIX 50 = Taenite 
0 0 


Mineral 


¢ The modulus values in the unit cell were computed based on two different 
published values of pyroxene. 


¢ The analysis shows significant difference in computed values of modulus. 


1. B. Kiefer, Elasticity of (Mg, Fe)SiO3 — Perovskite at High Pressures, Geophysical Research Letters, Volume 29, No. 11, 
2002. 

2. R. Liebermann, Elasticity Of Pyroxene-Garnet And Pyroxene-llmenite Phase Transformations In Germanates, Physics of the 
Earth and Planetary Interiors, 8 (1974) 361-374 
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Young’s Modulus for Stony Meteorites 


Comparison of Youngs Modulus 
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Youngs Modulus (GPa) 


These values were generated assuming the lower values for pyroxene. 
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Poisson’s Ratio for Stony Meteorites 


—— H-Type 
—— L-Type 
—— LL-Type 


Comparison of Poissons Ratio 
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Effect of Porosity 


Comparison of Youngs Modulus 


— —H Type 15% Matrix 

— —L Type 15% Matrix 

— — LL Type 15% Matrix 

— H Type 15% Matrix 10% Porosity 

——— L Type 15% Matrix 10% Porosity 
—— LL Type 15% Matrix 10% Porosity 
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The value of Young’s modulus go down significantly (~20%) with 10% 
porosity 
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Effect of Matrix 


Youngs Modulus Comparison for Mineral Constituent 10% and 15% Matrix 
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— 15% Matrix 
— 10% Matrix 
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5% change in matrix would cause ~15 GPa difference in modulus 
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Combine Effects of Different Pyroxene 
Modulus and Porosity 


H Type L Type and LL Type Comparison 


Lower pyroxene E 
10% porosity 


Higher pyroxene E, 
0% porosity 
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Validation of Meteorite Unit Model 


¢ Two materials are being tested as well as published data is 
referenced to validate the meteorite unit model 


— Flood Basalt from Pullman, Washington (dense rock with little 
porosity or cracks) 


— Tamdakht Meteorite (H-type ordinary chondrite, Fall in Morocco) 
¢ Following methods were used for validation 

— Uni-axial compression test 

— Acoustic Velocity and density measurements 

— Published data in literature 


SS 
PSK OCTOBER 23 - 27, 2016 
ase ose | SALT PALACE CONVENTION CENTER | SALT LAKE CITY, UTAH USA 


yr =~ -») dw 


Basalt - Unit Model Development 


BASALT MINERAL COMPOSITION % VOLUME 


Peak at 5% 
@ Plagioclase Ca-rich 


CT] Columbia River 


2 
@ Plagioclase Na-rict i meee 
@ Si-K20 Glass : Other Basalt 
@ Orthopyroxene § 

@ Pigeonite 


@ High-Ca Pyroxene 


@ Olivine 

@ Porosity co) 10 20 30 40 50 60 70 80 
Total Porosity ,in percent 

@ Matrix 


Figure 6. Histogram of total porosity of basalts (modified from Gu- 
zowski and others, 1982, E-2). 
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Unit Model Predictions for Basalt 


Columbia River Basalt Youngs Modulus 


— 10% Matrix 5% Porosity 
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Basalt Unit 
10% MATRIX & 5% POROSITY 
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Unit Model Validation - Basalt 


Comparison of Youngs Modulus 


= Published #1 Hi 


= Published #1 Low 

= Published #2 

= Unit Cell 10% Porosity 
= Unit Cell 5% Porosity 
m™ Average Acoustic Data 


Unit Cell values are close to published data. 


Youngs Modulus (GPa) 
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Tamdakht - Unit Model Development 


TAMDAKHT MINERAL COMPOSITION % PLAGIOCLASE Ca-Rich 
VOLUME 


PLAGIOCLASE Na-Rich 
@ Olivine 
@ Pyroxene 
@ Diopside 
@ Feldspar 
@ Triolite ORTHOPYROXENE 
@ Kamacite 
M@ Taenite 
mw Matrix 
@ Porosity 


Si-K20 Glass 


PIGEONITE 


A, 


High-Ca PYROXENE 


Olivine Tamdakht Unit: 15% 
Matrix & 10% Porosity 

MATRIX 

POROSITY 
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Unit Predictions: Tamdakht and Generic 
H-type Meteorite 


Comparison of Youngs Modulus 
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Tamdakht Measured Values — Acoustic Data 


Young's Modulus from Acoustic Data 


= Sample 1 - x axds 

= Sample 1 - y axis 
Sample 1 - y axis opposite 

= Sample 1 - z axis 

m Sample 2 - x axis 


= Sample 2 - y axis 
= Sample 2 - z axis 
a fl = Unit Cell 


Significantly lower values than unit cell model predictions (~100 
GPa) based on pure mineralogy. Next step is to introduce 
cracks/voids in the unit cell model and recomputed the values. 


Youngs Modulus (GPa) 


OCTOBER 23 — 27, 2016 
7 SALT PALACE CONVENTION CENTER | SALT LAKE CITY, UTAH USA 


Meteorite Unit Development: Morphology 


¢ While homogenization of different constituents help take care of the 
compositional variation, it does not address the influence on effective structural 
properties due to sharp crack and other morphological features in ground based 
meteorites. 


¢ 1000s of ordinary chondrites characterized as representative of their stony 
asteroid parent bodies. 


¢ The results indicate that, when present most of the fractures are (1) in the form 
of thin veins and usually show no obvious orientation, (2) in some cases veins 
radiated from a point of weakness, (3) occasionally veins have chicken-wire or 
a brick-wall network 
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As expected, presence of cavity greatly influenced the stress-strain relation and 


The stress-strain curve for models with and without cracks were analyzed. 
effective modulus. 


In order to investigate the influence of cracks 


homogenized meteorite properties. 


Strength Model Development - 
Compression tests 


¢ Uniaxial compression tests are being 
performed 

¢ Cylinders are chosen to prevent effects of 
sharp corners. 

¢ Each cylindrical sample will be strain- 
gaged to obtain accurate strain to failure, 
stress-strain curve as well as modulus 
values. 

¢ CT scans will be performed to provide 
morphological map. 

¢ The test will be simulated in FE model 

¢ Strength models will be developed based 
on results from compression tests 

¢ The data will be further used to develop 
Weibull parameters and scaling laws. 
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Preliminary Tests 


¢ 4samples (2 Basalt 1 cm cubes, 1 
Tamdakht 1 cm cube and 1 basalt 
rectangular sample were tested on 
Instron testing machine.) 


¢« Basalt shows classic brittle fracture. 
Samples broke at the sharp corners 


that led to lower strength values Basalt Tamdakht 
compared to published data as shown 
in the bar chart on right. Compressive Strength 


¢ Tamdakht sample didn't break at the 
sharp corners, instead we saw several a 
cracks developing at the surface and 
fracture was more gradual. 


¢ The data from one test shows lower 
strength compared to flood basalt. 50 
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Conclusions and Future Outlook 


We were successfully able to develop the meteorite unit models 
that account for complex mineralogy, porosity and matrix 
architecture. 


We were able to validate the meteorite unit model for terrestrial 
rock like flood Basalt. 


Cracks and voids will be introduced next, in the meteorite unit 
model based on the data from CT scans. 


Compression tests will be performed in near future to investigate 
stress-strain relationship, modulus, distribution of strength and 
strain to failure. 


The laboratory tests will be simulated to validate the meteorite 
unit models. 

After successful validation and verification of meteorite unit 
models the results will be used for material model inputs for full 
scale asteroid simulations. 

The results at different scales will also be used to develop Weibull 
parameters for predictions of break-up during entry. 
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